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Entanglement is one of the most intriguing features of quantum mechanics. It describes non-local correla¬ 
tions between quantum objects, and is at the heart of quantum information sciences. Entanglement is rapidly 
gaining prominence in diverse fields ranging from condensed matter to quantum gravity. Despite this generality, 
measuring entanglement remains challenging. This is especially true in systems of interacting delocalized parti¬ 
cles, for which a direct experimental measurement of spatial entanglement has been elusive. Here, we measure 
entanglement in such a system of itinerant particles using quantum interference of many-body twins. Lever¬ 
aging our single-site resolved control of ultra-cold bosonic atoms in optical lattices, we prepare and interfere 
two identical copies of a many-body state. This enables us to directly measure quantum purity, Renyi entangle¬ 
ment entropy, and mutual information. These experiments pave the way for using entanglement to characterize 
quantum phases and dynamics of strongly-correlated many-body systems. 


At the heart of quantum mechanics lies the principle of su¬ 
perposition: a quantum system can be in several states at the 
same time. Measurement on such a superposition state will 
exhibit randomness in the outcomes. This quantum random¬ 
ness is fundamental in nature, unlike classical randomness 
that arises when the observer has incomplete knowledge or 
ignores information about the system, as when throwing dice 
or flipping coins. In a many-body quantum system, quantum 
superposition between various possible configurations often 
results in a correlated randomness in the measurement out¬ 
comes of different parts of the system. These correlated sub¬ 
systems are then said to be entangled m. The non-local cor¬ 
relations between entangled subsystems prompted Einstein to 
describe entanglement as ‘spooky action at a distance’ 0, and 
were shown by Bell to be inconsistent with reasonable local 
theories of classical hidden variables 0. Later, it was real¬ 
ized that entanglement could be used as a resource to perform 
tasks not possible classically, with applications in computa¬ 
tion in 13, communication 0, and simulating the physics of 
strongly correlated quantum systems El. 

In few level quantum systems, entangled states have been 
investigated extensively for studying the foundations of quan¬ 
tum mechanics il and as a resource for quantum informa¬ 
tion applications 013. Recently, it was realized that the 
concept of entanglement has broad impact in many areas of 
quantum many-body physics, ranging from condensed mat¬ 
ter ITOl to high energy field theory ifTTl and quantum grav¬ 
ity ini. In this general context, entanglement is most often 
quantified by the entropy of entanglement m that arises in 
a subsystem when the information about the remaining sys¬ 
tem is ignored. This entanglement entropy exhibits qualita¬ 
tively different behavior than classical entropy and has been 
used in theoretical physics to probe various properties of the 
many-body system. In condensed matter physics, for exam¬ 
ple, the scaling behavior II3 of entanglement entropy allows 
distinguishing between phases that cannot be characterized by 
symmetry properties, such as topological states of matter Dl- 
[161 spin liquids iflTl [TSl . Entanglement entropy can be 
used to probe quantum criticality |[T9l and non-equilibrium 
dynamics EqIEI], and to determine whether efficient numer¬ 
ical techniques for computing many-body physics exist ||22l. 
Despite the growing importance of entanglement in theoreti- 
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Eigure 1. Bipartite entanglement and partial measurements. 

A generic pure quantum many-body state has quantum correlations 
(shown as arrows) between different parts. If the system is divided 
into two subsystems A and B, the subsystems will be bipartite en¬ 
tangled with each other when quantum correlations span between 
them (right column). Only with no bipartite entanglement present, 
the partitioned system I'i/'As) can be described as a product of sub¬ 
system states lip a) and |V’s) (left column). A path for measuring 
the bipartite entanglement emerges from the concept of partial mea¬ 
surements: ignoring all information about subsystem B (indicated as 
“Trace”) will put subsystem A into a statistical mixture, to a degree 
given by the amount of bipartite entanglement present. Finding ways 
of measuring the many-body quantum state purity of the system and 
comparing that of its subsystems would then enable measurements 
of entanglement. Eor an entangled state, the subsystems will have 
less purity than the full system. 


cal physics, current condensed matter experiments do not have 
a direct probe to detect and measure entanglement. Synthetic 
quantum systems such as cold atoms EH ED, photonic net¬ 
works ll25l . and some microscopic solid state devices ESI 
have unique advantages: their almost arbitrary control and 
detection of single particles, experimental access to relevant 
dynamical time scales, and isolation from the environment. 
In these systems, specific entangled states of few qubits, such 
as the highly entangled Greenberger-Horne-Zeilinger (GHZ) 
state El have been experimentally created and detected using 
witness operators ll2^ . However, entanglement witnesses are 
state specific. An exhaustive method to measure entanglement 
of an arbitrary state requires reconstructing the quantum state 
using tomography 1291 . This has been accomplished in small 






















2 


systems of photonic qubits lf30l and trapped ion spins OTl . but 
there is no known scheme to perform tomography for systems 
involving itinerant delocalized particles. With multiple copies 
of a system, however, one can use quantum many-body in¬ 
terference to quantify entanglement even in itinerant systems 

GUEIlEa. 

In this work, we take advantage of the precise control and 
readout afforded by our quantum gas microscope 0^ to pre¬ 
pare and interfere two identical copies of a four-site Bose- 
Hubbard system. This many-body quantum interference en¬ 
ables us to measure quantities that are not directly accessible 
in a single system, e.g. quadratic functions of the density ma¬ 
trix Ell m [331 |35}ill. Such non-linear functions can re¬ 
veal entanglement IH. In our system, we directly measure the 
quantum purity, Renyi entanglement entropy, and mutual in¬ 
formation to probe the entanglement in site occupation num¬ 
bers. 


BIPARTITE ENTANGLEMENT 

To detect entanglement in our system, we use a fundamen¬ 
tal property of entanglement between two subsystems (bipar¬ 
tite entanglement): ignoring information about one subsys¬ 
tem results in the other becoming a classical mixture of pure 
quantum states. This classical mixture in a density matrix 
p can be quantified by measuring the quantum purity, de¬ 
fined as Tr(p^). For a pure quantum state the density ma¬ 
trix is a projector and Tr(p^) = 1, whereas for a mixed 
state Tr(p^) < 1. In case of a product state, the subsys¬ 
tems A and B of a many-body system AB described by a 
wavefunction |'0 ab) (Fig- are individually pure as well, 
i.e. Tr(p^) = Tr(p|) = Tt:{p\b) = 1. Here the re¬ 
duced density matrix of A, pA = At:b{pab), where pab — 
\4’ab){'4'ab\ is the density matrix of the full system. Trg 
indicates tracing over or ignoring all information about the 
subsystem B. For an entangled state, the subsystems become 
less pure compared to the full system as the correlations be¬ 
tween A and B are ignored in the reduced density matrix, 
Tr(p^) = Tr(/9^) < Tr(p^g) = 1. Even if the many-body 
state is mixed (Tr(p^^) < 1), it is still possible to measure 
entanglement between the subsystems IT]. It is sufficient 1^ 
to prove this entanglement by showing that the subsystems are 
less pure than the full system, i.e. 

Tr(pA) < Tr(p^B), 

Tr(p|) < Av{p\b)- (1) 

These inequalities provide a powerful tool for detecting entan¬ 
glement in the presence of experimental imperfections. Fur¬ 
thermore, quantitative bounds on the entanglement present in 
a mixed many-body state can be obtained from these state pu¬ 
rities ll40l . 

Eq.Q can be framed in terms of entropic quantities ifTl I^ . 
A particularly useful and well studied quantity is the n-th or¬ 
der Renyi entropy, 

5„(A) = ^logTr(p:i). (2) 
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Figure 2. Measurement of quantum purity with many-body 
bosonic interference of quantum twins, a. When two N parti¬ 
cle bosonic systems that are in identical pure quantum states are in¬ 
terfered on a 50%-50% beam splitter, they always produce output 
states with even number of particles in each copy. This is due to the 
destructive interference of odd outcomes and represents a general¬ 
ized Hong-Ou-Mandel interference, in which two identical photons 
always appear in pairs after interfering on a beam splitter, b. If the 
input states pi and p 2 are not perfectly identical or not perfectly pure, 
the interference contrast is reduced. In this case the expectation value 
of the parity of particle number (Pi) in output i = 1, 2 measures the 
quantum state overlap between the two input states. For two identi¬ 
cal input states pi = p 2 , the average parity (Pi) therefore directly 
measures the quantum purity of the states. We only assume that the 
input states have no relative macroscopic phase relationship. 

From Eq. we see that the second-order (n = 2) Renyi 
entropy and purity are related by 52 (^) = —logTr(p^). 
S 2 ( 2 I) provides a lower bound for the von Neumann entangle¬ 
ment entropy 5yAr(A) = 5i(A) = —Tr(p^ logp^) exten¬ 
sively studied theoretically. The Renyi entropies are rapidly 
gaining importance in theoretical condensed matter physics, 
as they can be used to extract information about the “entan¬ 
glement spectrum” ED providing more complete informa¬ 
tion about the quantum state than just the von Neuman en¬ 
tropy. In terms of the second-order Renyi entropy, the suffi¬ 
cient conditions to demonstrate entanglement ID139) become 
52 (A) > 52(AB), and 52 (B) > 52(AB), i.e. the subsys¬ 
tems have more entropy than the full system. These entropic 
inequalities are more powerful in detecting certain entangled 
states than other inequalities like the Clauser-Home-Shimony- 
Holt (CHSH) inequality EH [33. 


MEASUREMENT OE QUANTUM PURITY 

The quantum purity and hence the second-order Renyi en¬ 
tropy can be directly measured by interfering two identical 
and independent copies of the quantum state on a 50%-50% 
beam splitter ETI [33[33[[Ml- For two identical copies of a 
bosonic Fock state, the output ports always have even parti- 
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cle number, as illustrated in Fig. |^. This is due to the de¬ 
structive interference of all odd outcomes. If the system is 
composed of multiple modes, such as internal spin states or 

(fc) 

various lattice sites, the total number parity Pi = YikPi 
is equal to unity in the output ports i = 1,2. Here the par¬ 
ity for mode k, = ±1 for even or odd number of parti¬ 
cles, respectively. The well known Hong-Ou-Mandel (HOM) 
interference of two identical single photons ll42ll is a special 
case of this scenario. Here a pair of indistinguishable pho¬ 
tons incident upon different input ports of a 50%-50% beam 
splitter undergoes bosonic interference such that both photons 
always exit from the same output port. In general, the average 
parity measured in the many-body bosonic interference on a 
beam splitter probes the quantum state overlap between the 
two copies {Pi) = Tr(/5i/52), where pi and p 2 are the density 
matrices of the two copies respectively and (...) denotes aver¬ 
aging over repeated experimental realizations or over identical 
systems, as shown in Fig. |^. Hence, for two identical sys¬ 
tems, i.e. for Pi = p 2 = p, the average parity for both output 
ports {i = 1,2) equals the quantum purity of the many-body 
state ET] [321133, 


(P,) = Tr(p2). (3) 

Equation Q represents the most important theoretical foun¬ 
dation behind this work - it connects a quantity depending on 
quantum coherences in the system to a simple observable in 
the number of particles. It holds even without fixed particle 
number, as long as there is no definite phase relationship be¬ 
tween the copies (Supplementary material). From Eqs. 0 
and Q, detecting entanglement in an experiment reduces to 
measuring the average particle number parity in the output 
ports of the multi-mode beam splitter. 

We probe entanglement formation in a system of interacting 
®^Rb atoms on a one dimensional optical lattice with a lattice 
constant of 680 nm. The dynamics of atoms in the lattice is 
described by the Bose-Hubbard Hamiltonian, 

= -1), (4) 

(i.i) * 


where a], Ui and Ui = a|ai are the bosonic creation and an¬ 
nihilation operators, and the number of atoms at site i, re¬ 
spectively. The atoms tunnel between neighboring lattice sites 
(indicated by (z, j)) with a rate J and experience an onsite re¬ 
pulsive interaction energy U. The Planck’s constant h is set 
to 1 and hence both J and U are expressed in Hz. The di¬ 
mensionless parameter (7/ J is controlled by the depth of the 
optical lattice. Additionally, we can superimpose an arbitrary 
optical potential with a resolution of a single lattice site by us¬ 
ing a spatial light modulator (SLM) as an amplitude hologram 
through a high resolution microscope (Supplementary mate¬ 
rial). This microscope also allows us to image the number 
parity of each lattice site independently li34l . 

To initialize two independent and identical copies of a state 
with fixed particle number N, we start with a low entropy 
2D Mott insulator with unity filling in the atomic limit ll34l 
and deterministically retain a plaquette of 2 x atoms while 
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Figure 3. Many-body interference to probe entanglement in op¬ 
tical lattices, a. A high resolution microscope is used to directly 
image the number parity of ultra cold bosonic atoms on each lat¬ 
tice site (raw images: green = odd, black = even). Two adjacent 
ID lattices are created by combining an optical lattice and potentials 
created by a spatial light modulator (SLM). We initialize two iden¬ 
tical many-body states by filling the potentials from a low entropy 
2D Mott insulator. The tunneling rates Jx, Jy can be tuned indepen¬ 
dently by changing the depth of the potential, b. The atomic beam 
splitter operation is realized in a tunnel coupled double well poten¬ 
tial. An atom, initially localized in one of the wells, delocalizes with 
equal probability into both the wells by this beam splitter. Here, we 
show the atomic analog of the HOM interference of two states. The 
joint probability "P)!, 1) measures the probability of coincidence de¬ 
tection of the atoms in separate wells as a function of normalized 
tunnel time Jyt, with the single particle tunneling Jy = 193(4) Hz. 
At the beam splitter duration {Jyt — 1 /8) bosonic interference leads 
to a nearly vanishing P{1,1) corresponding to an even parity in the 
output states. This can be interpreted as a measurement of the purity 
of the initial Fock state, here measured to be « 0.90(4). The data 
shown here are averaged over two independent double wells. The 
blue curve is a maximum likelihood fit to the data, and the error- 
bars reflect 1 a statistical error, c. When two copies of a product 
state, such as the Mott insulator in the atomic limit are interfered on 
the beam splitter, the output states contain even number of particles 
globally (full system) as well as locally (subsystem), indicating pure 
states in both. d. On the other hand, for two copies of an entangled 
state, such as a superfluid state, the output states contain even num¬ 
ber of particles globally (pure state) but a mixture of odd and even 
outcomes locally (mixed state). This directly demonstrates entangle¬ 
ment. 
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removing all others (Supplementary material). This is illus¬ 
trated in Fig. [^. The plaquette of 2 x iV atoms contains two 
copies (along the y-direction) of an iV-atom one-dimensional 
system (along the ^-direction), with = 4 in this figure. The 
desired quantum state is prepared by manipulating the depth 
of the optical lattice along x, varying the parameter U / Jx 
where Jx is the tunneling rate along x. A box potential cre¬ 
ated by the SLM is superimposed onto this optical lattice to 
constrain the dynamics to the sites within each copy. During 
the state preparation, a deep lattice barrier separates the two 
copies and makes them independent of each other. 

The beam splitter operation required for the many-body in¬ 
terference is realized in a double well potential along y. The 
dynamics of atoms in the double well is likewise described 
by the Bose-Hubbard hamiltonian, Eq. 0- A single atom, 
initially localized in one well, undergoes coherent Rabi oscil¬ 
lation between the wells with a Rabi frequency of J = Jy (os¬ 
cillation frequency in the amplitude). At discrete times during 
this evolution, t = tss = the atom is delocal¬ 

ized equally over the two wells with a fixed phase relation¬ 
ship. Each of these times realizes a beam splitter operation, 
for which the same two wells serve as the input ports at time 
t = 0 and output ports at time t = Ibs- Two indistinguish¬ 
able atoms with negligible interaction strength {U/Jy 1) 
in this double well will undergo interference as they tunnel. 
The dynamics of two atoms in the double well is demon¬ 
strated in Eig. in terms of the joint probability V{1, 1 ) 
of finding them in separate wells versus the normalized time 
Jyt. The joint probability V{1, 1 ) oscillates at a frequency 
of 772(16) Hz = 4Jy, with a contrast of about 95(3)%. At 
the beam splitter times, 7^(1,1) « 0. The first beam split¬ 
ter time, tJy = I is used for all the following experiments, 
with 7^(1,1) « 0.05(2). This is a signature of bosonic inter¬ 
ference of two indistinguishable particles ||4^ l44l akin to the 
photonic HOM interference ll42ll . This high interference con¬ 
trast indicates the near-perfect suppression of classical noise 
and fluctuations and includes an expected 0.6% reduction due 
to finite interaction strength (U/ Jy « 0.3). The results from 
this interference can be interpreted as a measurement of the 
quantum purity of the initial Eock state as measured from the 
average parity (Eq.0), {Pi) = 1 — 2 x 7^(1,1) = 0.90(4), 
where i = 1, 2 are the two copies. 


ENTANGLEMENT IN THE GROUND STATE 

The Bose-Hubbard model provides an interesting system to 
investigate entanglement. In optical lattice systems, a lower 
bound of the spatial entanglement has been previously esti¬ 
mated from time-of-flight measurements 1461 and entangle¬ 
ment dynamics in spin degrees-of-freedom has been inves¬ 
tigated with partial state reconstruction El. Here, we di¬ 
rectly measure entanglement in real space occupational par¬ 
ticle number in a site-resolved way. In the strongly interact¬ 
ing, atomic limit of U/ Jx ^ 1, the ground state is a Mott 
insulator corresponding to a Eock state of one atom at each 
lattice site. The quantum state has no spatial entanglement 
with respect to any partitioning in this phase - it is in a prod¬ 


uct state of the Eock states. As the interaction strength is re¬ 
duced adiabatically, atoms begin to tunnel across the lattice 
sites, and ultimately the Mott insulator melts into a superfluid 
with a fixed atom number. The delocalization of atoms create 
entanglement between spatial subsystems. This entanglement 
originates HSjjSOll from correlated fluctuations in the number 
of particles between the subsystems due to the super-selection 
rule that the total particle number in the full system is fixed, 
as well as coherence between various configurations without 
any such fluctuation. 

To probe the emergence of entanglement, we first prepare 
the ground state of Eq. Q in both the copies by adiabatically 
lowering the optical lattice potential along x. Then we freeze 
the tunneling along x without destroying the coherence in the 
many-body state and apply the beam splitter along y. Einally, 
we rapidly turn on a very deep 2D lattice to suppress all tun¬ 
neling and detect the atom number parity (even = 1, odd = -1) 
at each site. We construct the parity of a spatial region by mul¬ 
tiplying the parities of all the sites within that region. The av¬ 
erage parity over repeated realizations measures the quantum 
purity, both globally and locally according to Eq. 0^ enabling 
us to determine the second-order Renyi entropy globally and 
for all possible subsystems. In the atomic Mott limit (Eig. |^), 
the state is separable. Hence, the interference signal between 
two copies should show even parity in all subsystems, indi¬ 
cating a pure state with zero entanglement entropy. Towards 
the superfluid regime (Eig. |^), the buildup of entanglement 
leads to mixed states in subsystems, corresponding to a finite 
entanglement entropy. Hence, the measurement outcomes do 
not have a pre-determined parity. Remarkably, the outcomes 
should still retain even global parity, indicating a pure global 
state. Higher entropy in the subsystems than the global sys¬ 
tem cannot be explained classically and demonstrates bipartite 
entanglement. 

Experimentally, we find exactly this behavior for our two 
4-site Bose-Hubbard systems (Eig. 4). We observe the emer¬ 
gence of spatial entanglement as the initial atomic Mott insu¬ 
lator melts into a superfluid. The measured quantum purity 
of the full system is about 0.6 across the Mott to superfluid 
crossover, corresponding to a Renyi entropy, S 2 {AB) « 0.5. 
The measured purity deep in the superfluid phase is slightly 
reduced, likely due to the reduced beam splitter fidelity in 
presence of increased single sites occupation number, and any 
residual heating. The nearly constant global purity indicates 
high level of quantum coherence during the crossover. Eor 
lower interaction strength U j Jx (superfluid regime) we ob¬ 
serve that the subsystem Renyi entropy is higher than the full 
system, S' 2 (A) > S 2 {AB). This demonstrates the presence of 
spatial entanglement in the superfluid state. In the Mott insu¬ 
lator regime (U / Jx 3> 1), 5*2(A) is lower than S 2 {AB) and 
proportional to the subsystem size, consistent with a product 
state. 

In these experiments, we post-select outcomes of the ex¬ 
periment for which the total number of atoms detected in both 
copies is even. This constitutes about 60% of all the data, 
and excludes realizations with preparation errors, atom loss 
during the sequence, or detection errors (Supplementary ma¬ 
terial). The measured purity is consistent with an imperfect 
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Figure 4. Entanglement in the ground state of the Bose-Hubhard model. We study the Mott insulator to superfluid transition with four 
atoms on four lattice sites in the ground state of the Bose-Hubbard model, Eq. (|^. a. As the interaction strength U/ Jx is adiabatically reduced 
the purity of the subsystem A (green and blue, inset), Tr(p^), become less than that of the full system (red). This demonstrates entanglement in 
the superfluid phase, generated by coherent tunneling of bosons across lattice sites. In terms of the second-order Renyi entanglement entropy, 
S 2 {A) = — log Tr(p^), the full system has less entropy than its subsystems in this state. In the Mott insulator phase {Uj Jx 1) the full 

system has more Renyi entropy (and less purity) than the subsystems, due to the lack of sufficient entanglement and a contribution of classical 
entropy. The circles are data and the solid lines are theory calculated from exact diagonalization. The only free parameter is an added offset, 
assumed extensive in system size and consistent with the average measured entropy in the full system, b. Second-order Renyi entropy of all 
possible bi-partitioning of the system. For small U/Jx, all subsystems (data points connected by green and blue lines) have more entropy 
than the full system (red circles), indicating full multipartite entanglement between the four lattice sites. The residual entropy in the Mott 
insulating regime is from classical entropy in the experiment, and extensive in the subsystem size. Right: The values of all Renyi entropies of 
the particular case of U/Jx ~ 1 are plotted, to demonstrate spatial multipartite entanglement in this superfluid. 


beam splitter operation alone, suggesting significantly higher 
purity for the many-body state. The measured entropy is thus 
a sum of an extensive classical entropy due to the imperfec¬ 
tions of the beam splitter and any entanglement entropy. 

Our site resolved measurement simultaneously provides in¬ 
formation about all possible spatial partitionings of the sys¬ 
tem. Comparing the purity of all subsystems with that of the 
full system enables us to determine whether a quantum state 
has genuine spatial multipartite entanglement where every site 
is entangled with each other. Experimentally we find that this 
is indeed the case for small U/Jx (Fig. |^). In the super¬ 
fluid phase, all possible subsystems have more entropy than 
the full system, demonstrating full spatial multipartite entan¬ 
glement between all four sites ll^ |45]| . In the Mott phase 
(U/Jx ^ 1), the measured entropy is dominated by extensive 
classical entropy, showing a lack of entanglement. 

By measuring the second-order Renyi entropy we can cal¬ 
culate other useful quantities, such as the associated mutual 
information Iab = 5'2(A) + 82 ( 8 ) — S 2 {AB). Mutual in¬ 
formation exhibits interesting scaling properties with respect 
to the subsystem size, which can be key to studying area laws 
in interacting quantum systems ED. In some cases, such as 
in the ‘data hiding states ’ EH, mutual information is more 
appropriate than the more conventional two point correlators 
which might take arbitrarily small values in presence of strong 
correlations. Mutual information is also immune to present 
extensive classical entropy in the experiments, and hence is 
practically useful to experimentally study larger systems. In 
our experiments (Fig. we find that for the Mott insula¬ 


tor state (U/Jx ^ 1), the entropy of the full system is the 
sum of the entropies for the subsystems. The mutual informa¬ 
tion Iab ~ 0 for this state, consistent with a product state in 
the presence of extensive classical entropy. At U/ Jx ~ 10, 
correlations between the subsystems begin to grow as the sys¬ 
tem adiabatically melts into a superfluid, resulting in non-zero 
mutual information, Iab > 0. 

It is instructive to investigate the scaling of Renyi entropy 
and mutual information with subsystem size iniEii since 
in larger systems they can characterize quantum phases, for 
example by measuring the central charge of the underlying 
quantum field theory ifTTI . Figure shows these quantities 
versus the subsystem size for various partitioning schemes 
with a single boundary. For the atomic Mott insulator the 
Renyi entropy increases linearly with the subsystem size and 
the mutual information is zero, consistent with both a product 
state and classical entropy being uncorrelated between vari¬ 
ous sites. In the superfluid state the measured Renyi entropy 
curves are asymmetric and first increase with the system size, 
then fall again as the subsystem size approaches that of the full 
system. This represents the combination of entanglement en¬ 
tropy and the linear classical entropy. This non-monotonicity 
is a signature of the entanglement entropy, as the entropy for a 
pure state must vanish when the subsystem size is zero or the 
full system. The asymmetry due to classical entropy is absent 
in the mutual information. 

The mutual information between two subsystems comes 
from the correlations across their separating boundary. For 
a four site system, the boundary area ranges from one to three 
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Figure 5. Renyi Mutual information in the ground state. Contribution from the extensive classical entropy in our measured Renyi entropy 
can be factored out by constructing the mutual information Iab = S 2 {A) + S 2 {B) — S 2 {AB). Mutual information takes into account all 
correlations (HI between the subsystems A and B. a. We plot the summed entropy S 2 {A) + S 2 {B) (in blue, green and light blue corresponding 
to the partitions shown) and the entropy of the full system S 2 {AB) (in red) separately. Mutual information is the difference between the two, as 
shown by the arrow for a partitioning scheme. In the Mott insulator phase {U(Jx ^ 1) the sites are not correlated, and Iab ~ 0. Correlations 
start to build up for smaller U/Jx, resulting in a non-zero mutual information. The theory curves are from exact diagonalization, with added 
offsets consistent with the extensive entropy in the Mott insulator phase, h. Classical and entanglement entropies follow qualitatively different 
scaling laws in a many-body system. Top - In the Mott insulator phase classical entropy dominates and S 2 {A) and 52(B) follow a volume 
law- entropy increases with the size of the subsystem. The mutual information Iab ~ 0. Bottom - Sa, Sb show non-monotonic behavior, 
due to the dominance of entanglement over classical entropy, which makes the curves asymmetric. Iab restores the symmetry by removing 
the classical uncorrelated noise, c. Top - More correlations are affected (red arrow) with increasing boundary area, leading to a growth of 
mutual information between subsystems. The data points are for various partitioning schemes shown in Fig. Bottom- Iab as a function 
of the distance d between the subsystems shows the onset and spread of correlations in space, as the Mott insulator adiabatically melts into a 
superfluid. 


for various partitioning schemes. Among those schemes with 
a single boundary maximum mutual information in the super¬ 
fluid is obtained when the boundary divides the system sym¬ 
metrically (Fig. Bi). Increasing the boundary size increases 
the mutual information, as more correlations are interrupted 
by the partitioning (Fig.[^). 


Mutual information also elucidates the onset of correlations 
between various sites as the few-body system crosses over 
from a Mott insulator to a superfluid phase. In the Mott in¬ 
sulator phase (U/Jx ^ 1) the mutual information between 
all sites vanish (Fig. [^, bottom). As the particles start to 
tunnel only the nearest neighbor correlations start to build up 
{U/Jx ~ 12) and the long range correlations remain negli¬ 
gible. Further into the superfluid phase, the correlations ex¬ 
tend beyond the nearest neighbor and become long range for 
smaller U/Jx- These results suggest disparate spatial behav¬ 
ior of the mutual information in the ground state of an uncor¬ 
related (Mott insulator) and a strongly correlated phase (su¬ 
perfluid). For larger systems this can be exploited to identify 
quantum phases and the onset of quantum phase transitions. 


NON-EQUILIBRIUM ENTANGLEMENT DYNAMICS 

Away from the ground state, the non-equilibrium dynam¬ 
ics of a quantum many-body system is often theoretically in¬ 
tractable. This is due to the growth of entanglement beyond 
the access of numerical techniques such as the time depen¬ 
dent Density Matrix Renormalization Group (DMRG) theory 
1531 m. Experimental investigation of entanglement may 
shed valuable light onto non-equilibrium quantum dynamics. 
Towards this goal, we study a simple system: two particles 
oscillating in a double well ll43l |55]| . This non-equilibrium 
dynamics are described by be Bose-Hubbard model. The 
quantum state of the system oscillates between unentangled 
(particles localized in separate wells) and entangled states in 
the Hilbert space spanned by |1,1), |2,0) and |0, 2). Here, 
|m,n) denotes a state with m and n atoms in the two sub¬ 
systems (wells), respectively. Starting from the product state 
11,1) the system evolves through the maximally entangled 
states |2,0) -F |0, 2) ± |1,1) and the symmetric HOM-like 
state |2, 0) -F |0, 2). In the maximally entangled states the 
subsystems are completely mixed, with a probability of 1 /3 
to have zero, one, or two particles. The system then returns 
to the initial product state |1,1) before re-entangling. In our 
experiment, we start with a Mott insulating state (U/ Jx ^ 1), 
and suddenly quench the interaction parameter to a low value, 
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Figure 6. Entanglement dynamics in quench Entanglement dy¬ 
namics of two atoms in two sites after a sudden quench of the Hamil¬ 
tonian from a large value of U/ Jx to U/ Jx ~ 0.3, with Jx ~ 210 
Hz. Here, ‘quench time’ refers to the duration that the atoms spend 
in the shallow double well, after the initial sudden quench. The sys¬ 
tem oscillates between Mott insulator like state (I) and quenched su¬ 
perfluid states (II, III). The growth of bipartite entanglement in the 
superfluid state is seen by comparing the measured Renyi entropy of 
the single site subsystem (blue circles) to that of the two site full sys¬ 
tem (red circles). The solid lines are the theory curves with vertical 
offsets to include the classical entropy introduced by experimental 
imperfections. 


Uj Jx ~ 0.3. The non-equilibrium dynamics is demonstrated 
(Fig. 1^ by the oscillation in the second-order Renyi entropy 
of the subsystem, while the full system assumes a constant 
value originating from classical entropy. This experiment also 
demonstrates entanglement in HOM-like interference of two 
massive particles. 


SUMMARY AND OUTLOOK 

In this work, we perform a direct measurement of quan¬ 
tum purity, the second-order Renyi entanglement entropy, and 
mutual information in a Bose-Hubbard system. Our measure¬ 
ment scheme does not rely on full density matrix reconstruc¬ 
tion or the use of specialized witness operators to detect en¬ 
tanglement. Instead, by preparing and interfering two iden¬ 
tical copies of a many-body quantum state, we probe entan¬ 


glement with the measurement of only a single operator. Our 
experiments represent an important demonstration of the use¬ 
fulness of the many-body interference for the measurement of 
entanglement. It is straight forward to extend the scheme to 
fermionic systems ll5^ and systems with internal degrees of 
freedom ll^ . By generalizing the interference to n copies of 
the quantum state ll35l . arbitrary observables written as as n- 
th order polynomial function of the density matrix, e.g. n > 2 
order Renyi entropies, can be measured. 

With modest technical upgrades to suppress classical fluc¬ 
tuations and residual interactions, it should be possible to fur¬ 
ther improve the beam splitter hdelity enabling us to work 
with signihcantly larger systems. Mutual information is an 
ideal tool for exploring these larger systems as it is insensi¬ 
tive to any residual extensive classical entropy. For high en¬ 
tropy of a sub-system, corresponding to low state purity, the 
number of measurements required to reach a desired preci¬ 
sion is high. However, in contrast to tomographic methods, 
this scheme would not require additional operations for larger 
systems. Moreover, the single site resolution of the micro¬ 
scope allows us to simultaneously obtain information about 
all possible subsystems, to probe multipartite entanglement. 

For non-equilibrium systems, entanglement entropy can 
grow in time (indefinitely in inhnite systems). This leads to in¬ 
teresting many-body physics, such as thermalization in closed 
quantum systems 122]. The long time growth of entanglement 
entropy is considered to be a key signature of many-body lo¬ 
calized states IMI arising in presence of disorder. The ability 
to measure the quantum purity for these systems would allow 
experimental distinction of quantum fluctuations and classical 
statistical fluctuations. 

More generally, by starting with two different quantum 
states in the two copies this scheme can be applied to measure 
the quantum state overlap between them. This would provide 
valuable information about the underlying quantum state. For 
example, the many-body ground state is very sensitive to per¬ 
turbations near a quantum critical point. Hence, the overlap 
between two ground states with slightly different parameters 
(such as ( 7 / J in the Bose-Hubbard hamiltonian) could be used 
as a sensitive probe of quantum criticality ll58] . Similarly the 
overlap of two copies undergoing non-equilibrium evolution 
under different perturbations can be used to probe temporal 
correlation functions in non-equilibrium quantum dynamics. 
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SUPPLEMENTARY MATERIAL 

I. MEASURING ENTANGLEMENT ENTROPY WITH 
QUANTUM INTEREERENCE 

The quantification of entanglement requires the measure¬ 
ment of non-linear functionals of a quantum state p, such as 


the n-th order Renyi entropy = — In Tr(p”) HI. A general 
scheme to measure Tr(p") is to measure the shift operator 
acting on n-copies of the many-body system. The shift oper¬ 
ator Vn re-orders the quantum states when acting on a collec¬ 
tion of n states, 

K|t/’l)|t/’2) • • ■ \tpn) = |V'n)l'*/'l) • ■ • IV'n-l). (S.l) 

It can be shown that Tr(p”) = Tr(I4p®”) El. 

We focus on the experimentally relevant case of n = 2. 

The shift operator is then simply the SWAP operator which 
exchanges any two quantum states: 

V' 2 (|'!/'i) ® IV' 2 )) = IV' 2 ) ® IV'l) (S-2) 

Two successive applications of the SWAP operator leave the 
system unchanged, V 2 = 1. Therefore V 2 has eigenvalues 
±1, corresponding to subspaces of the 2-copy system that 
are symmetric or antisymmetric with respect to the state ex¬ 
change. The SWAP operator may act on individual modes 
(e.g. lattice sites) or the entire quantum system, and op¬ 
erations on different modes commute. The following short 
proofEl 

Tr (V 2 P 1 0 P 2 ) = Tr ^¥3 ^ Pm |i) 01 ® |k) (1| j 

= |k) 01 |i) (1| j 

= H pf^Pki^kjSa = = Tr (P1P2) 

ijkl ik 

(5.3) 

shows that the overlap of two quantum states pi and p 2 is 
given by the expectation value of the SWAP operator on the 
product space of the two states. Consider from now on the 
case where the two state are identical {pi = P 2 = p), then 
the expectation value of V 2 gives the purity Tr(p^). Further 
if we have two copies of a pure state then Tr(p^) = 1 = 
Tr(V2/9® p), hence the combined 2-copy state is symmet¬ 
ric and can be expressed in the symmetric basis comprised 
of states 

|vac) > with p,q = 0,1,2 ,... 

i,3 ) 

(5.4) 

where 0 ^ 2 ) i '^he creation operation of mode i in copy 1(2). 

If the two copies undergo a discrete Fourier transformation of 
the form (for simplicity dropping the mode indices) 

(a| + oY)!\/2 — y a\ 

{a\- a\)/s/2^ a\ (S.5) 

then the basis states in Eq |S.4| will end up having 2 x pi parti¬ 
cles in mode i of copy 2. In other words a symmetric state, as 
is the case for 2 pure identical copies, will always have even 


10 


arbitrary pure state 


copy I 1^) 


l«i> 


no phase relation (number state) 




jni = even) Pi 


copy II |W) 


|u2 = even) l^'v) 


|n2 = even) P2 


(-P2) = {V2) = Trp^ ^ {Pi) 


{Pi) = {Vi, = = (Pi) 


Figure 7. Beamsplitter for many-body interference. Left: With 
the beamsplitter operation and proper phase shift operations, one can 
directly measure quantum purity by measuring the average parity in 
output port 2 of the beamsplitter. For pure identical incident states, 
the atom number is always even in output 2. Right: In the exper¬ 
iment, we interfere states with well-defined particle number N or 
subsystems of such states. No macroscopic phase relationship exists 
between the input states, and the phase shifts in the input/output ports 
have no physical significance. Both outputs are equivalent and may 
be used to measure the expectation value (V 2 ) of the SWAP operator. 


number of particles in copy 2 after the transformation. The 
symmetric and anti-symmetric subspaces of the SWAP opera¬ 
tor are identified by the parity of atom number in copy 2 after 
a discrete Fourier transform, and the average parity directly 
measures the state purity, (P 2 ) = Trp^. 

Our microscope experiments then allows us to probe en¬ 
tanglement in an optical lattice by comparing the local purity 
Tr/9^(5)^ to the global purity Trp^ for a system partitioned 
into subsystems A and B. The entanglement is quantified by 
the entropy of the reduced density matrix pA = Tyb{p), and 
the measured purity Trp^ gives directly the 2nd order Renyi 
entropy S 2 = — InTrp^. This scheme is proposed in |j2l and 
made explicit for measurements with beamsplitter operations 
in optical lattices in El and El. giving Renyi entropy of 
arbitrary order n. 


Using controlled tunneling in a double-well potential, we 
can implement the beamsplitter transformation for bosonic 
atoms (see next section): 


I — y (fl| + io,2')I 

{ia\ + a\)lV2 


4 


(S.6) 


where a 7r/2-phase (f) is associated with each tunneling event 
across the double-well. Note that this transformation is not 


equivalent to the Fourier transform in Eq. (S.5 1 . It’s easy to 


verify that the Fourier transform is realized with the following 
protocol sequence of the beamsplitter operation and relative 
phase shift operations: 

1. A —7r/2 phase shift (—i) on copy 2 


2. The beamsplitter operation in Eq. (S. 61 

3. Another —7r/2 phase shift on copy 2 


The inclusion of the additional phase shifts are important 
to correctly map the symmetric (antisymmetric) eigenstates 


of the SWAP operator onto states with even (odd) atom num¬ 
ber parity in output port 2 of the beamsplitter. The resulting 
protocol is valid for measuring purity of any general many- 
body state. In the classical limit where the incident states 
are two identical coherent states with well-defined identical 
phases, the inclusion of the proper phase factors in input 2 
ensures that the states interfere destructively in output 2. In 
this port, the total number of bosons is always zero and there¬ 
fore even, so the measured parity {V) = 1 correctly gives 
Pi = P 2 and Trp^ = 1. This situation is analogous to the 
interference of two phase-stabilized laser beams on a 50/50 
beamsplitter, which may result in zero intensity in one out¬ 
put for the correct choice of incident phases. Away from the 
classical limit, for example as the input states become number 
squeezed states with decreasing uncertainty in atom number 
but increasing fluctuation in phase, atoms start to appear in 
output port 2 after the protocol but only in parrs (even parity) 
as long as the input states remain pure and identical. 

The protocol also works when there is no global phase re¬ 
lationship between the interfering many-body states. Such as 
in our current experiments when the two copies are prepared 
each as an independent quantum state with a fixed number of 
atoms, so there is no well-defined phase. There are also no 
defined phases when the incident states to the beamsplitter 
are subsystems partitioned out of bigger systems. In either 
case, step 1 of the above protocol has no physical meaning 
in the absence of a defined phase and might be omitted from 
the experiment without changing the resulting state after 
the transformation. The in-situ fluorescence imaging of our 
microscope detects parity of the atom number on each lattice 
site which is phase-insensitive, so step 3 is also redundant. 
The beamsplitter operation in the double-well alone is thus 
sufficient to implement the mapping of SWAP operator eigen¬ 
states onto states with even or odd atom number parity. The 
two output ports are then equivalent and the purity measure¬ 
ments may be obtained from the atom numbers on either side 
of the double-well after the many-body interference sequence. 


II. IMPLEMENTATION OF THE BEAMSPLITTER 
OPERATION 


A. Projected double-well potentials 


In addition to a square lattice, optical potentials are gener¬ 
ated by projecting light onto the atoms using a digital micro¬ 
mirror device (DMD). The DMD is used as an amplitude holo¬ 
gram in a Fourier plane of our high resolution imaging system 
so that wavefronts with arbitrary phase and amplitude may be 
created with single-site resolution 0. We use blue-detuned 
coherent light at A = 762 nm to generate a potential with a 
double-well profile along x and a smoothed flat top profile 
along y. 
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X (lattice site) 

Figure 8. Double-well potential for the beamsplitter. The inten¬ 
sity profile of the projected potential (Eq |S.7[ gray short-dashed), the 
lattice (gray long-dashed), and the combined potential for the beam 
splitter operation (black solid). Also shown are sketches of the am¬ 
plitude of the ground band Wannier wavefunctions (blue, purple) in 
each well at the beamsplitter depth. 



Figure 9. Rabi oscillations in the double-well. A single particle 
is initialized on the right side of the doublewell and oscillates coher¬ 
ently between the two wells with fitted tunneling rate J = 193(4) FIz 
and contrast of 95(1)%. 


_ (x-1.5)^ _ _(x + 1.6)2 

V[x,y) = Vdw[e 0 . 95 ^ — 0.52e ^+e 0 . 95 ^ 


X arctan 


y + 18 
5.5 


— arctan 


t/-18 

5.5 


(S.7) 


where x and y are in units of the lattice spacing and Vdw is 
the potential depth of the projected double-well. 

The beamsplitter operation is realized by controlled tunnel¬ 
ing in the combined potential of the above projected poten¬ 
tial and a shallow at-lattice, as depicted in Fig.[^ We choose 
depths Vdw = 1-7 and Viatt = 2 E^, for which we observe 
tunneling rate J = 193(4) Hz during the beam-splitter oper¬ 
ation (Fig.|^, in reasonable agreement with a band structure 
calculation predicting J = 170 Hz. The discrepancy is likely 
due to uncertainly in the lattice depth, which is calibrated 
using amplitude modulation spectroscopy at Viatt ~ iOEr- 
Here En = 1240 Hz is the recoil energy of the optical lattice. 
In the beamsplitter potential, the energy gap to the first excited 
band is » 1.3 kHz, and states outside the ground band do not 
contribute significantly to the dynamics. 


B. Sources of error for the beamsplitter 

a. Potential imperfections The leading order imperfec¬ 
tion of the projected double-well potential are imperfect zero- 
crossings in the electric field, resulting in energy offsets be¬ 
tween the two sides of the double-well. At the double-well 
depth for our beamsplitter operation, we observe offsets of 
50 Hz or less, which do not significantly affect the Hong-Ou- 
Mandel (HOM) interference contrast (see Fig.[T0|i. 



U/J 


A=0 J 
A=0.6 J 
A=0.8 J 


Figure 10. Fidelities of the beamsplitter operation. Finite inter¬ 
actions and energy offsets due to imperfections in the double-well 
potential reduce the Hong-Ou-Mandel interference contrast, as mea¬ 
sured by the probability to detect even atom numbers at the beam¬ 
splitter time tJ = J. For a beamsplitter operation starting with 
one atom on each side of the double-well and typical experimental 
parameters J = 240 Hz, U = 70 Hz and offset A = 50 Hz (corre¬ 
sponding to ( 7 / J = 0.3 and A/J — 0.2 ), the interference contrast 
is expected to be reduced by ~ 0.6%. This calculation does not take 
the effects of higher bands into account. 


b. Alignment stability The successful loading of atoms 
from the lattice into the double-well potential is sensitive to 
long-term and shot-to-shot position drifts between the lattice 
and the double-well. We minimize such drifts by imaging the 
lattice and double-well potential at the end of each experimen¬ 
tal run and feeding back on the position of the double-well 
with a piezo-actuated mirror. We achieve a relative position 
stability of 0.04 sites RMS or less. To lowest order the posi¬ 
tion drift creates an energy offset between the two sides of the 
combined double-well potential. At the chosen depths for the 
beamsplitter operation, a relative shift of 0.1 sites leads to an 
offset of « 20Hz. 

c. Interaction effects Interactions during the beamsplit¬ 
ter operation potentially reduce the HOM interference con- 
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Mott insulator Number state preparation Many-body evoiution Beam-splitter Fluor, 

preparation (2xN-site) (2indep. IxN) operation imaging 



Figure 11. Experimental sequence. Schematic showing the ramps 
of the X- and y- lattices and the projected potential from the DMD. 
The profiles of the DMD potentials are sketched in the direction of 
interest, while the other direction always has a smoothed flattop pro¬ 
file aross the region of interest. Ramps are exponential in depth as a 
function of time. See text for details. 


trust. We minimize interactions by performing all experiments 
in a weak out-of-plane confinement of = 27r x 800 Hz. 
During the beamsplitter operation we achieve an interaction 
of [7 = 27r X 70 Hz (measured with photon-assisted tunnel¬ 
ing in a deep double-well and extrapolated to lower depths), 
corresponding to u = U/J = 0.3. This residual interac¬ 
tion reduces the HOM interference contrast by ~ 0.6% (see 

Fig.[T0|. 

d. Coherent admixture of higher bands Interactions of 
two particles on the same site distort the particles’ wavefunc- 
tions and coherently admix higher bands. This wavefunc- 
tion is thus different from that of a single particle, restor¬ 
ing some distinguishability to the bosonic atoms. The dom¬ 
inant contribution of higher bands occurs in the z-direction, 
along which the confinement is weakest, and the second ex¬ 
cited band is admixed to the wavefunction. The admixture is 
e « iirf = ( 2 x 2 ;J 8 Tok )^ = 0-2%- HOM interference 
contrast is thus reduced by less than 1%. 


III. EXPERIMENTAL SEQUENCE 

Our experiments start with a single-layer two-dimensional 
Mott insulator of ®^Rb atoms in a deep lattice (14 — = 

A5Er) with 680 nm spacing as described in previous work. 
The following sequence is illustrated in FigfTT] 

State preparation We deterministically prepare a plaque- 
tte of 2 X 2 or 2 X 4 atoms from a Mott insulator with unity oc¬ 
cupancy. We first superimpose onto the deep lattice an optical 
potential with a double-well profile of depth « bOEr along y 
and a smooth flattop profile along x, and subsequently turn off 
the lattice in the y-direction. The two troughs of the double¬ 
well are aligned with two neighboring lattice sites so only two 


rows of atoms are trapped, while all other atoms leave the sys¬ 
tem along tubes in the y-direction. A blue-detuned Gaussian 
beam with waist w « 50 ym and A = 755 nm provides the 
anti-confinement to remove atoms outside the double-well ef¬ 
ficiently in 40 ms. The y-lattice is then ramped back to its 
initial depth and the double-well removed, leaving a block of 
width 2 sites and length ~ 10 sites populated with one atom 
each. The above procedure is then repeated with a double- or 
quadruple-well potential along the x direction, leaving a deter¬ 
ministically loaded block of 2 x 2 or 2 x 4 atoms in the lattice. 
The lattices and double-well potentials are ramped smoothly 
to avoid heating the atoms to excited bands of the lattice. 

At the end of the state preparation sequence, the fidelity 
of unit occupancy is 94.6(2)% per site, limited primarily by 
the fidelity of the initial Mott insulator and losses during the 
state preparation. We verify independently that defects are 
predominantly empty, not doubly occupied sites. 

Evolution in independent copies For studying the ground 
state entanglement using the 2x4 block (Figure. 4 & 5 in the 
main text), we turn on an optical potential with two narrow 
Gaussian peaks separated by four lattice sites along the x di¬ 
rection and flat-top along y. This confines the atoms inside the 
4-site ’’box potential”. The x lattice is then ramped down adi- 
abatically to various final depths from 26 to 377^. The ramp in 
depth is exponential in time with a duration of 200 ms from 45 
to iEr- The y-lattice is kept at AbEr so that tunneling along 
y is negligible and the two copies evolve independently. 

For quench dynamics using the 2x2 block, we use a 
double-well potential along x with Vdw = 277^ to prevent 
atoms from leaving the 2-site system. The x lattice is ramped 
from 4577r to 277^ in 0.75 ms and held for a variable time. 
The presence of the double-well slightly modifies the values 
of U and J from values in a lattice only, yielding U/ J = 0.3 
during the hold time. 

Beamsplitter operation and imaging Right before the 
beamsplitter operation, the a;-lattice is ramped back to 4577,. 
in 0.75 ms to suppress tunneling within each copy. A double¬ 
well potential along y is superimposed onto the lattice. The 
y-lattice is then ramped down to 277^ in 0.5 ms and atoms are 
allowed to tunnel in independent double-wells between the 
two copies for 0.34 ms, implementing the beam splitter trans¬ 
formation before the y-lattice is returned to its initial depth of 
AbEr in 0.5 ms. 

Subsequently, we pin the atoms in a deep lattice and obtain 
fluorescence images with single-lattice-site resolution. Our 
detection scheme is sensitive only to the parity of the site oc¬ 
cupation number. 


IV. DATA ANALYSIS 

Post-selection Before data analysis we post-select out¬ 
comes of the experiment for which the total number of 
atoms detected within the plaquette (2 x 2 or 2 x 4 sites) is 
even. Outcomes outside this subset of data indicate either 
state preparation errors, atom loss during the sequence, or 
detection errors. We furthermore reject all realizations for 
which we detect atoms in the 20 x 20 block surrounding 
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Postselect total even No postselection 




Figure 12. Renyi entropy of the 4-site system and its subsystems with and withont postselection. The postselection process removes 
classical entropy and reduces the entropy of the full system (red) from ~ 0.9 to ~ 0.5. Note that even without postselection the entropy of 
the half-chain (green) rises above the full system entropy, indicating bipartite entanglement. Theory curves are exact diagonalizations shifted 
vertically by the mean classical entropy per site calculated from the full system entropy. 


the region of interest, most likely corresponding to atoms 
being lost from the plaquette during the sequence. Note that a 
combination of multiple errors (e.g. failure to load atoms on 
two sites) may lead to an unsuccessful run of the experiment 
being included in the post-selected data. 

For experiments on the 2x2 plaquette we typically retain 
80% of the data, and 60% for the 2 x 4 plaquette. 

Calculating Purity and Entropy For the full many-body 
state and each subsystem of interest we calculate Podd, the 
probability of measuring a total odd number of atoms after 
the beamsplitter operation within the post-selected data. The 
quantum mechanical purity and second-order Renyi entropy 
are then given by 

Tr(p^) = 1 - 2podd (S.8) 

52(p) = -log(Tr(p2)) (S.9) 

We average the calculated purity Tr(p^) over both copies 
and over equivalent partitions. For instance, the single-site en¬ 
tropy reported in Fig. 4a of the main text is the mean over the 
first and last site of each copy of the 4-site system. From the 
variance of the parity in each subsystem and the covariance 
between subsystems we obtain the statistical standard error 
of the mean parity, taking into account possible correlations 
between regions. The reported error bars are the quadrature 
sum of the statistical error and the standard deviation of mean 
parities over the averaged regions. This procedure accounts 
for residual inhomogeneities between the copies and along the 
chains. 

Errorbars mU / J correspond to a typical uncertainty in the 
optical lattice depth of ±2%. 


Full system purity For the 2x4 plaquette, the initial 
state purity is reduced from 1 due to the presence of thermal 
holes in the initial Mott insulating state. Assuming all even 
sites are holes, the loading statistics for the 2x4 plaquette are: 


N atoms 

loading probability p{N) 

8 

0.66(1) 

7 

0.27(1) 

6 

0.052(4) 


i.e. the postselected subset of total even data contains 
n nro°^n ge = 7% of cvcnts with 6 atoms total. The inclu- 

0 . 052 + 0.66 

sion of outcomes with 6 atoms reduces the purity of the initial 
state to 0.94, corresponding to a Renyi entropy of 0.06. The 
expected full system purity in the Mott insulator state is thus 
limited by the finite 95(3)% fidelity of the beamsplitter op¬ 
eration on each site and approximately given by the product 
of individual purities, P = ~ 0.90"^ « 0.66, in good 

agreement with the experimentally measured purity in Fig. 4a. 

Fitting procedure To determine the contrast of single¬ 
particle Rabi oscillations (Fig|^ and HOM-interference 
(Fig. 3b in main text) we use a Bayesian inference for the 
fit to the measured parity, which is more robust than a least- 
squares fit in situations where error probabilities are small and 
the visibility close to 1. This approach prevents unphysical 
fits that extend past the possible bounds of the measurement 
and appropriately adjusts the error bars for points measured 
to lie near the physical bound. For each time point, we re¬ 
port the median and the lcr(68%) confidence interval of a 
Beta-distribution f3{p,m,N) for m successful outcomes in 
N experimental runs. The fitted sine curves in Fig. 1 are 
maximum-likelihood fits to the Beta-distributions at each time 
point, which are determined by maximizing the product of all 
the Beta-distributions where the fitted curve samples them. El 
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